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We investigate the Andreev transport through double quantum dot Cooper pair splitters with 
ferromagnetic leads. The analysis is performed with the aid of the real-time diagrammatic technique 
in the sequential tunneling regime. We study the dependence of the Andreev current, the differential 
conductance and the tunnel magnetoresistance on various parameters of the model in both the linear 
and nonlinear response regimes. In particular, we analyze the spin-resolved transport in the crossed 
Andreev reflection regime, where a blockade of the current occurs due to enhanced occupation of the 
triplet state. We show that in the triplet blockade finite intradot correlations can lead to considerable 
leakage current due to direct Andreev reflection processes. Furthermore, we find additional regimes 
of current suppression resulting from enhanced occupation of singlet states, which decreases the 
rate of crossed Andreev reflection. We also study how the splitting of Andreev bound states, 
triggered by either dot level detuning, finite hopping between the dots or magnetic field, affects 
the Andreev current. While in the first two cases the number of Andreev bound states is doubled, 
whereas transport properties are qualitatively similar, in the case of finite magnetic field further 
level splitting occurs, leading to a nontrivial behavior of spin-resolved transport characteristics, and 
especially that of tunneling magnetoresistance. Finally, we discuss the entanglement fidelity between 
split Cooper pair electrons and show that by tuning the device parameters hdelity can reach unity. 

PACS numbers: 73.23.-b,73.21.La,74.45.+c,72.25.-b 


I. INTRODUCTION 

Quantum dots coupled to superconducting and normal 
leads provide very promising systems to study the inter¬ 
play between the superconducting correlations and the 
mesoscopic electronic transportSuch hybrid nanos¬ 
tructures have recently attracted a lot of attention due to 
the possibility to control and split the Cooper pairs 
When the quantum dot is attached to normal and su¬ 
perconducting lead and for voltages smaller than the su¬ 
perconducting energy gap A, the current flows through 
the system due to Andreev reflection.— More specifically, 
transport occurs then through sub-gap, Andreev bound 
states (ABS), which were recently probed experimentally 
using bias spectroscopy^^”— For three-terminal systems, 
e.g. with one superconducting and two normal leads, the 
Cooper pair, when leaving the superconductor, can en¬ 
ter either to the same normal lead or can be split when 
the two electrons forming Cooper pair end in different 
leads. The former process is known as direct Andreev 
reflection (DAR), whereas the latter one is referred to 
as crossed Andreev reflection (CAR). Usually both pro¬ 
cesses contribute to the Andreev current, however, under 
certain conditions, by properly changing device param¬ 
eters, one can tune the contributions due to CAR and 
DAR processes or even suppress one of them. This can 
be obtained in the case when the normal leads are fer¬ 
romagnetic. Then, in the antiparallel magnetic configu¬ 
ration of the device, for leads with large degree of spin 
polarization, only CAR processes contribute, since each 
lead supports electrons of opposite spin.— 

Another interesting system, which allows for con¬ 
trollable manipulation of Cooper pairs, can be made 
of double quantum dots (DQDs). Recent experiments 


have shown that DQDs can work as Cooper pair beam 
splitters, whose operation can be controlled by gate 
voltages A"— In contrast to single quantum dot hybrid 
systems, the DQD setup allows to study pure CAR trans¬ 
port regime by considering suitable system’s parameters. 
Specifically, in real double quantum dots the intradot 
Coulomb repulsion can be much larger than other energy 
scales,— which for a wide range of applied bias volt¬ 
ages prevents double occupancy of each dot. As a con¬ 
sequence, the direct Andreev reflection processes, which 
require simultaneous transfer of two electrons with oppo¬ 
site spins by the same dot, become suppressed. However, 
Andreev reflection processes can still occur through CAR 
processes, in which the electrons forming the Cooper pair 
are transferred simultaneously through the two dots. An¬ 
other advantage of DQDs is the possibility of independent 
level tuning of each individual dot. Due to this ability, it 
has been shown experimentally that Cooper pair splitting 
can be dominant on resonance, whereas out of resonance 
elastic cotunneling processes dominate^^i^ 

Transport properties of DQD Cooper pair beam 
splitters have already been addressed in several 
publications^”— which, among others, addressed 
the problem of coherence and entanglement of split 
Cooper pairs and their probing^”— as well as the 
noise correlations^^*^ and Cooper pair microwave 
spectroscopy ] These investigations were performed 
for DQD Cooper pair beam splitters with nonmagnetic 
leads. However, because using ferromagnetic leads can 
be important to estimate entanglement between split 
electrons,— providing comprehensive study of transport 
properties of DQD Cooper pair splitters with ferromag¬ 
netic contacts seems desirable. The analysis of Andreev 
transport through such systems is thus the goal of the 
present paper. Furthermore, DQDs coupled to a su- 
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FIG. 1. (Color online) Schematic of a double quantum 
dot Cooper pair splitter with ferromagnetic leads. The left 
(right) dot is coupled to the left (right) lead with the cou¬ 
pling strength FJ (F^) and each dot is coupled to a common 
superconducting lead with coupling Ff and F;^ for the left and 
right dot, respectively. The magnetizations of the ferromag¬ 
netic leads can form either parallel or antiparallel magnetic 
configuration. The level energy and Coulomb repulsion in 
dot i are denoted by Si and Ui, while t and Ulr describe the 
hopping and the Coulomb correlations between the two dots. 


merit fidelity between split electrons forming Cooper pair 
and show that, depending on parameters, fidelity can 
reach unity. 

The paper is organized in the following way: Sec. 
II contains the description of the DQD model and the 
method used in calculations. The numerical results and 
their discussion in the crossed Andreev reflection regime 
are presented in Sec. III. In Sec. IV we analyze how the 
splitting of Andreev bound states affects the transport 
properties. The next section is devoted to the analysis of 
transport in the full parameter space where both CAR 
and DAR processes are present. The behavior of entan¬ 
glement fidelity on bias and gate voltages is studied in 
Sec. VI and, finally, the paper is concluded in Sec. VII. 


II. THEORETICAL FRAMEWORK 

A. Model Hamiltonian 


per conductor and to two ferromagnetic contacts, can 
exhibit a considerable tunnel magnetoresistance (TMR) 
and generate large spin current. Therefore, such nanos¬ 
tructures are also interesting for spin nanoelectronics 
and understanding their magnetoresistive properties is 
of great importance. Transport properties of quan¬ 
tum dots attached to ferromagnetic leads have already 
been broadly investigated both experimentally^iri^ and 
theoretically!^-— However, spin-resolved transport prop¬ 
erties of hybrid dots, consisting of quantum dots cou¬ 
pled to ferromagnetic and superconducting leads, have 
been mainly studied in the case of single quantum 
dots^-— while the case of double quantum dots is 
largely unexplored 

In this paper, we therefore investigate spin-dependent 
Andreev transport through two single-level quantum 
dots, coupled to one superconducting and two ferromag¬ 
netic leads. Our analysis is performed with the aid of the 
real-time diagrammatic technique in the lowest-order ex¬ 
pansion with respect to the coupling to ferromagnetic 
leads, while the coupling to superconductor can be ar¬ 
bitrarily strong. First, by assuming infinite correlations 
in the dots, we analyze the pure CAR regime where the 
triplet blockade of the current occurs !2^ We then show 
that even very large but finite intradot correlations can 
lead to considerable leakage current in the triplet block¬ 
ade. We thoroughly study the behavior of the Andreev 
current, differential conductance and TMR, deriving ap¬ 
proximate zero-temperature formulas for TMR in appro¬ 
priate transport regimes. We also analyze the effect of fi¬ 
nite splitting of Andreev bound states on transport prop¬ 
erties, discussing the splitting caused by either level de¬ 
tuning, finite hopping between the dots or finite magnetic 
field. For finite correlations in the dots, we study trans¬ 
port properties in the full parameter space and identify 
additional transport regimes where the current suppres¬ 
sion occurs. At the end, we also consider the entangle- 


We consider double quantum dot Cooper pair split¬ 
ter, which is schematically displayed in Fig. [TJ It con¬ 
sists of two single-level quantum dots, each attached to 
its own ferromagnetic lead, and both coupled to a com¬ 
mon s-wave superconductor. The magnetizations of fer¬ 
romagnetic leads are assumed to form either parallel (P) 
or antiparallel (AP) configuration. Switching between 
these two configurations can be obtained upon applying 
a small external magnetic field Bg. We assume that this 
field is so small that it does not lead to the splitting of the 
dot’s level, neither affects it the superconducting phase. 
The total system is modeled by the following effective 
Hamiltonian: 

H= Hp + Hs + Hdqd + Ht, ( 1 ) 

^=L,R 

where the first term, describes the left {/3 = 

L) and right (/3 = R) ferromagnetic electrodes in 
the noninteracting quasiparticle approximation, = 
Ek<T ek/3<TCk/3<T'^k/3<T- Here, 4^^ (ck/3a) is the creation (an- 
nihilation) operator of an electron with the wave vector 
k and spin a in the lead /3, whereas 5k/3cr denotes the 
corresponding single-particle energy. The second term in 
Eq. ^ describes the s-wave BCS superconducting (S) 
lead in the mean field approximation 

Hs = y^gk5o-Ck5o-<^k5cr + A (CkS'tC-kS't + h.C.) (2) 

kcr k 

with e\^sa denoting the relevant single-particle energy 
and A standing for the order parameter of the super¬ 
conductor, which is assumed to be real and positive. 
The third term of the Hamiltonian m describes the 
two single-level quantum dots and acquires the follow- 




3 


ing form: 

Hdq]j — ^ ^ I ^ ^ -f- B^Siz “h UiTii^Tii^ I 

i=L,R \ cr / 

+ ^ ULRriLanRa' + t '^{d\^^dRa- + d^^^dRa), ( 3 ) 

a,a' cr 

where d\^ creates a spin-cr electron in dot i of energy 5^, 
^icr = d\^dicj^ Ul (Ur) is the Coulomb correlation energy 
of the left (right) dot, and Bz denotes external magnetic 
field in units of gjHB = ^ with Siz = — nii)/2. Urr 

and t stand for the interdot Coulomb repulsion and the 
hopping between the dots, respectively. 

The last term of the Hamiltonian describes tunneling 
of electrons between the leads (L, i?, S) and the two dots 

= E E (4) 

kcr i=L,R 

+EE (^ikcr<^kS'cr^*cr + h.C.), 

kcr i=L,R 

with (^iko-)’ i = denoting the relevant tun¬ 

neling matrix elements between the two dots and ferro¬ 
magnetic leads (the superconducting lead). In the fol¬ 
lowing, we assume that these matrix elements are k and 
a independent, and The cou¬ 

pling of the dots to respective ferromagnetic leads can be 
parametrized by, T^ = 27r|H*ppf, where pf is the spin- 
dependent density of states of lead i. Within the wide 
band approximation these couplings become energy inde¬ 
pendent and constant. Introducing the spin polarization 
of the z-th lead, Pi = {pt - Pi )/{pt + Pi where p+ {p~) 
is the spin majority (minority) density of states, the cou¬ 
plings can be written in the form, T^ = (1 -\-(jpi)Ti, with 
T^ = (rj -\-V\)/2. Generally, each dot can be coupled 
to its lead with different strength and the two leads can 
have different spin polarizations, here, however, we re¬ 
strict our analysis to symmetric systems and note that 
the presented results are also qualitatively valid for sys¬ 
tems with weak asymmetry in the couplings. We thus 
assume Pl = Pr = P and Tr = Tr = T/2. More¬ 
over, we also assume that the dots’ levels are degenerate 
El = £r = £ and the dots’ Coulomb energies are equal, 
Ul = Ur = I/, unless stated otherwise. 

Since in this paper we are only interested in Andreev 
transport, we can take the limit of an infinite supercon¬ 
ducting gap, A ^ oo. Then, the quantum dot system 
coupled to the superconducting lead can be described by 
the effective Hamiltonian^ 

= Hdqd - E W {4Al + h-c-) (5) 

i=L,R 

where Tf^ = yT^T^- The superconducting proxim¬ 
ity effects are included in the last two terms of Eq. (|5|). 


The first term describes local proximity effects on each 
dot and arises due to direct Andreev reflection, whereas 
the second term describes creation of nonlocal entangled 
states between the two dots. These nonlocal correlations 
are responsible for crossed Andreev reflection. The effec¬ 
tive pair potential Vf {i = L, R) is the coupling strength 
between the i-th dot and superconducting electrode and 
acquires the form, Tf = 27r\Vf\‘^ps^ where ps denotes 
the density of states of the superconductor in the normal 
state. We assume that the couplings between the dots 
and superconductor are equal, Tf = T^ = T^. 

The device is biased in the following way: The elec¬ 
trochemical potential of the superconducting lead is as¬ 
sumed to be grounded, ps = 0, see Fig. [TJ while the 
potentials of the left and right leads are kept the same, 
Pl = Pr = P = eV. In this way the net current between 
the left and right ferromagnetic lead is zero. In the follow¬ 
ing, we use the convention that for positive bias, eV > 0, 
the Cooper pairs tunnel to the superconductor, while for 
negative bias, eE < 0, the Cooper pairs are extracted 
from the superconducting electrode. 

We would like to note that while the assumption of 
infinite superconducting gap allows us to exclude normal 
tunneling processes and study only the Andreev trans¬ 
port, it needs to be taken with some care. This is be¬ 
cause in real systems the gap can be large, but is clearly 
finiteHowever, for relatively low bias voltages, as 
considered in this paper, one can expect normal tunnel¬ 
ing processes to be negligible and, thus, the assumption 
of large superconducting energy gap is reasonable. 

B. Method 

In order to calculate the transport characteristics of 
the considered system, we employ the real-time diagram¬ 
matic technique (RTDT)^“— which is based on system¬ 
atic perturbation expansion of the reduced density ma¬ 
trix and operators of interest with respect to the coupling 
strength T. Within the RTDT, in the stationary state the 
reduced density matrix p can be found fro m^^i^^ 

E^x,x'^x'=0, (6) 

x' 

where the elements of fho self-energy matrix W 

describe transitions between the states |x) and |x') on 
the Keldysh contour, with |x) denoting the many-body 
eigenstate of = ^xIx)- P is the vector 

of diagonal density matrix elements = (x|p|x), which 
can be found from Eq. (j6j) together with the normaliza¬ 
tion condition. The current flowing from the ferromag¬ 
netic lead i can be found fro m^^i^^ 

7i = |lV{W^'P}, (7) 

where W^* denotes the modified self-energy matrix W, 
which takes into account the number of electrons trans¬ 
ferred through the junction i. 
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To find the occupation probabilities and the cur¬ 
rent, we perform the perturbation expansion of the self¬ 
energies, occupation probabilities and the current with 
respect to the coupling strength to ferromagnetic leads T. 
In our studies, we consider the weak coupling regime and 
take into account only the first-order tunneling processes, 
which correspond to sequential tunneling. We first deter¬ 
mine the self-energies using the respective diagrammatic 
rules^i^ and then calculate the occupation probabilities 
and the current by using Eqs. and ©■ 

We also note that when T ^ T^, as considered in the 
present paper, and taking into account only the lowest- 
order tunneling processes in the coupling to ferromag¬ 
netic leads, the reduced density matrix becomes diagonal 
in the eigenbasis of the effective HamiltonianThis is 
why Eq. includes only the diagonal elements of the re¬ 
duced density matrix. Moreover, we would like to empha¬ 
size that while the perturbation expansion with respect 
to the coupling strength T is performed, no assumption 
on the strength of the Coulomb correlation parameters is 
imposed, and they are treated in an exact way. Eurther- 
more, the assumption of the weak coupling regime implies 
that the Kondo temperature of the system is exponen¬ 
tially small. Thus, at temperatures considered in calcu¬ 
lations, the correlations leading to the Kondo effect are 
irrelevant and do not need to be taken into account 


C. Quantities of interest 

The main quantity of interest is the current flowing 
through the system due to Andreev reflection processes. 
By calculating the currents /l and Ir flowing through the 
left and right junctions, the total current flowing into the 
superconductor can be simply obtained from the Kirch- 
hoff’s law 

Is = II Ir, ( 8 ) 

together with the corresponding differential conductance, 

Gs = dls/dV. 

Since the normal leads are ferromagnetic, the Andreev 
current depends on the magnetic configuration of the de¬ 
vice, which is assumed to be either parallel or antiparal¬ 
lel. We thus also calculate the tunnel magnetoresistance 
associated with the change of magnetic configuration of 
the system, which is defined asi^ 

tAP _ tP 

TMR = \p ^ , (9) 

where and denote the Andreev current flowing 
into the superconductor in the parallel and antiparal¬ 
lel magnetic configurations, respectively. Note that this 
definition is opposite to that in the case of the Julliere 
modeli^ This is because for hybrid quantum dots with 
superconducting and ferromagnetic leads, the Andreev 
current in the antiparallel configuration is usually larger 
than that in the parallel configuration^i^ 


In the following we present and discuss the results on 
the Andreev transport through DQD Cooper pair split¬ 
ters with ferromagnetic leads obtained within the sequen¬ 
tial tunneling approximation. We systematically study 
the behavior of the Andreev current, the associated dif¬ 
ferential conductance and the TMR in both the linear 
and nonlinear response regimes, exploring basically the 
whole parameter space of the considered model. In par¬ 
ticular, to analyze the transport regime where CAR pro¬ 
cesses are dominant, we first assume infinite Coulomb 
correlations on the dots, so that DAR processes are to¬ 
tally suppressed. We then relax this condition and allow 
for finite correlations in the dots to study the transport 
properties basically in the whole parameter space. We 
also analyze the effects of finite hopping between the dots, 
nonzero detuning of DQD levels and finite external mag¬ 
netic held. In addition, we calculate the entanglement 
fidelity between the split electrons forming Cooper pair. 

III. RESULTS IN THE CROSSED ANDREEV 
REFLECTION REGIME 

When the charging energy of each dot is much larger 
than Coulomb correlations between the dots and the ap¬ 
plied bias voltage, the rate of direct Andreev reflection 
is suppressed and only CAR processes are possible. This 
condition is in fact one of the main requirements for the 
system to work as Cooper pair beam splitterTherefore 
to analyze the Andreev transport in the case when DAR 
processes are absent we now assume U ^ oo. For the 
sake of simplicity of the following discussion, let us also 
at this point neglect the hopping term assuming t = 0 
(the role of finite hopping will be considered further). 

In the limit of infinite intradot correlations, double oc¬ 
cupation of each dot is forbidden and only nine states 
(out of 16) of the effective Hamiltonian (|5j) are relevant 
for transport. These are the following states: empty 
state |0,0), four singly occupied states |cr, 0), |0,cr), and 
four doubly occupied states |cr, cr') for cr, cr' where 

\a,/3) = \a)L\l3)R^ with |0)i and |cr)^ denoting empty 
and singly occupied states of dot i. Due to the hop¬ 
ping and additional terms in Eq. (j5j) resulting from the 
proximity effect, the Hamiltonian is not diagonal in the 
above basis. Diagonalizing the effective Hamiltonian, one 
finds a new basis consisting of the following states: four 
singly occupied states |cr, 0), |0,cr), three triplet states 
\To) = (l4-,t) + \ \To) = IcTjCr), and two states 

being linear combinations of empty state and singlet 
state, \S) = (|i,t) “ We note that the 

triplet states become decoupled from the superconduc¬ 
tor since Cooper pairs consist of two electrons with com¬ 
pensated spin, i.e., Spair = 0, and there is no coupling 
between the singlet and triplet states. The correspond¬ 
ing eigenenergies of the eigenstates given by Eq. (pR|) 
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FIG. 2. (Color online) The absolute value of the Andreev 
current calculated for (a) the parallel {!§) and (b) antiparallel 
(Is^) magnetic configurations as a function of detuning 6 = 
^s + Ulr and the applied bias eV. The parameters are: Ts = 
0.5, T = 0.015, r = 0.01, t = 0, = 0, with Ulr = 1 

the energy unit, and p = 0.5. The current is plotted in units 
lo = eV/h. 


are E± = (5/2 ± 5^, where S = 2e Ulr denotes de¬ 
tuning betwee n the sing let states and the empty state, 
while 2eA = \/(5^ + 2r| measures the energy difference 
between the states |+) and |—). 

The Andreev bound state (ABS) energies are defined 
as^ 

EtT = + 2r|, (11) 

where (a,/3 = ±. These energies are the excitation ener¬ 
gies between doublet and singlet states of the double dot 
decoupled from the ferromagnetic leads. 

Due to the assumption, ^ oo, double occupancy 
of each dot is forbidden and direct Andreev tunneling 
becomes totally suppressed. The only way to transfer 
charge between the double dot and the superconductor is 
by crossed Andreev reflection, the process which involves 
two electrons with opposite spins coming from different 


ferromagnetic leads. For eV > 0^ the two electrons tun¬ 
nel to the superconductor, while for eV < 0, the Cooper 
pairs are extracted from the superconducting electrode 
and entangled pairs of electrons are transmitted to ferro¬ 
magnetic leads (each electron ends in different lead). 


A. Andreev current and differential conductance 

In Fig. [2] we show the dependence of the absolute value 
of the Andreev current on the applied bias eV and the 
detuning parameter 6 = 2s E Ulr for the parallel and 
antiparallel magnetic configurations of the system. At 
low bias, the Andreev current is generally suppressed 
due to th e Coulomb blockade, except for two values of 
6, 1^1 = — 2F|, where = E^^^ = 0 and the 

corresponding Andreev bound states are at resonance. 
With increasing the bias, the current starts flowing once 
\eV\ > As the Andreev reflection becomes 

optimized for parameters corresponding to particle-hole 
symmetry point, the Andreev current reaches maximum 
value for small detuning, ^ ^ 0. These two features are 
similar to those observed in a tree-terminal system in¬ 
cluding a single quantum dot^^ However, in the present 
case the Andreev current exhibits a striking difference, 
as it does not reveal the symmetry with respect to the 
bias reversal, which has been present for the single dot 
system. In the present case the absolute value of the 
Andreev current reveals a very strong asymmetry with 
respect to the sign change of the bias. As can be seen in 
Fig. [21 for positive bias, for which the double dot becomes 
occupied by two electrons, the Andreev current ceases to 
flow. In fact, for eV > |(^ + Ulr)/2\, the double dot is 
in the triplet state, which explains the vanishing of the 
Andreev current, as the symmetry of the triplet state 
does not match the symmetry of the s-wave supercon¬ 
ductor. The blockade region is therefore independent of 
the magnetic configuration of the system - the Andreev 
current stops flowing due to the triplet blockade in both 
the parallel and antiparallel alignments, see Fig. |2l 

Figure [3] shows the dependence of the Andreev differ¬ 
ential conductance Gs = dis /dV on the bias voltage eV 
and detuning parameter 6. The sudden drop of the An¬ 
dreev current around the bias voltage eV ^ {6 ^ Ulr)/ 2 
for 6 /Ulr > —1 leads to the appearance of a pronounced 
negative differential conductance, see Fig. |3l which is 
present in both magnetic configurations. On the other 
hand, outside the triplet blockade the differential conduc¬ 
tance reveals positive peaks whenever the electrochemi¬ 
cal potential of ferromagnetic leads crosses one of the 
Andreev levels. 

Except for the asymmetry due to the triplet block¬ 
ade, one can also notice another asymmetry of both the 
current and differential conductance with respect to the 
bias reversal, which is visible for large detuning S, see e.g. 
S/Ulr = 2 in Figs. |2]and|3l This asymmetry especially 
reveals in the differential conductance when comparing 
the intensity of the low-bias peaks, i.e. the peak associ- 
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FIG. 3. (Color online) The differential conductance Gs — 
dIs/dV of the Andreev current in the parallel (Gs) and an- 
tiparallel (Gs^) magnetic configurations as a function of de¬ 
tuning 6 and applied bias voltage eV. Parameters are the 
same as in Fig. [2l 


ated with the Andreev level for eV > 0 (and (5 > 0) 

with the amplitude of the maximum associated with the 
level for eV < 0 (and (5 > 0), see Fig.[3l To under¬ 

stand this effect let us consider the differential conduc¬ 
tance for large value of detuning parameter djU^R = 2. 
The Andreev processes can occur in the system if the 
occupation of states |+) and/or |—) is finite. For large 
detuning the state |+) is high in energy and does not play 
any role in the considered bias voltage regime. The state 
which is relevant for Andreev transport is the state |—). 
One should also note that for large value of S the state 
|—) contains only relatively small admixture of the singlet 
state, cf. Eq. m, while the empty state is mostly occu¬ 
pied. This generally leads to small values of the current 
for large 6. 

For bias voltages such that E^^^ < eV < there 

are no Andreev levels in the transport window and the 
current flowing into/out of the superconductor is sup¬ 
pressed. When eV crosses one of those levels, either for 
positive or negative bias, the current starts to flow and a 



eV/U,, 


FIG. 4. (Golor online) The Andreev current (a) in the parallel 
(solid line) and antiparallel (dashed line) magnetic configura¬ 
tion and the corresponding differential conductance (b) as a 
function of the bias voltage calculated for detuning parameter 
SIUlr = 2. The other parameters are the same as in Fig. [2l 


peak appears in the differential conductance, see Figs. [2] 
and [3l It can be seen that the system becomes more 
transparent for eV < 0, when the Cooper pairs are ex¬ 
tracted from the superconductor, than for eE > 0, when 
one injects electron pairs into superconducting lead. This 
is because for positive bias voltage singly occupied states 
become populated and, since the singlet state is required 
for the current to flow into the superconductor, the An¬ 
dreev current is decreased. More specifically, when pass¬ 
ing the energy 5 (note that 5 = UlrI‘^ for S = 2Ulr) the 
probability of finding the double dot in state |—) becomes 
strongly suppressed at the cost of enhanced occupation of 
one-electron states, decreasing the Andreev current. On 
the other hand, for negative bias voltage the singly occu¬ 
pied states play little role in transport and the current is 
then larger compared to the case of eE >0. 

For a better visualization of the behavior of the An¬ 
dreev current and related differential conductance in dif¬ 
ferent transport regimes we also plot the cross-sections 
of transport characteristics for two values of detuning 
parameter S. Figure [U displays the bias dependence of 
the current and differential conductance in both mag¬ 
netic configurations calculated for 6/Ulr = 2. The 
current as a function of eV exhibits well-defined steps 
corresponding to consecutive Andreev bound states be¬ 
ing active in transport, while the differential conduc¬ 
tance shows the respective peaks. One can also see that 
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the differential conductance exhibits negative value due 
to the triplet blockade, which occurs for bias voltage 
eV ^ {S-\-Ulr)/‘2 and for 6 /Ulr > —1. Interestingly, the 
differential conductance in the parallel magnetic align¬ 
ment reveals additional negative differential conductance, 
which is not related to the triplet blockade, see Figs.[3l^a) 
andllKb). This negative differential conductance develops 
for S > — 2r| and for bias voltages eV ^ . 

The effect of negative differential conductance in the 
parallel configuration can be explained bearing in mind 
that formation of Cooper pairs involves two electrons 
with opposite spins. In the parallel alignment there are 
more electrons with one spin orientation than with the 
other one, thus, the rate of electron pairs is determined 
by the density of states of minority carriers. When the 
double dot starts to be occupied by odd number of elec¬ 
trons (here, singly occupied), the occupation probability 
of electrons with spin-up orientation increases whereas 
that of electrons with spin-down decreases. With fur¬ 
ther increase of the bias voltage the occupation of the 
spin-up level becomes greatly enhanced, whereas that of 
spin-down level becomes strongly suppressed, giving rise 
to nonequilibrium spin accumulation. As a consequence 
of spin accumulation, the Andreev current becomes also 
suppressed, which reveals as negative values in the differ¬ 
ential conductance. For reversed bias voltages the spin 
accumulation becomes irrelevant since the double dot is 
in the state |—), being rather empty with small admix¬ 
ture of singlet state. Thus, Cooper pairs can be more 
easily extracted from the superconductor compared to 
the opposite bias polarization. 

The Andreev current and differential conductance as 
a function of bias voltage in the absence of detuning are 
shown in Fig. [5l Now, one can clearly see that while for 
negative bias the current displays typical steps accompa¬ 
nied with peaks in dlsjdV^ for positive bias voltage the 
current first increases but then drops and becomes fully 
suppressed due to the triplet blockade, see Fig. ISfa). The 
associated negative differential conductance is clearly vis¬ 
ible in both magnetic configurations, see Fig. [5fb). 


B. Tunnel magnetoresistance 



FIG. 5. (Color online) The same as in Fig. |4] calculated for 
detuning parameter (5 = 0. 


TMR 



5/Ulr 


To observe more subtle differences between the parallel 
and antiparallel magnetic configurations one needs to use 
quantity which is more sensitive to a change of magnetic 
alignment of ferromagnetic leads. In Fig. [6] we present 
the dependence of TMR on detuning S and the bias 
voltage eV. In the region determined by the equation 
^ l(^ + ^Li?)/2|, where the current ceases to flow due 
to the triplet blockade, the TMR becomes indeterminate. 
This region is marked by white area in Fig. [6l Moreover, 
the TMR is strongly suppressed for 
and for the bias voltage Ef^^ < eV < E^^^. For this 
transport regime however the first-order processes are 
suppressed and to obtain correct value of TMR higher- 
order tunneling events should be considered. As has been 


FIG. 6. (Color online) The tunnel magnetoresistance TMR as 
a function of detuning S and the bias voltage eV. The white 
region indicates the range of parameters where the TMR is 
undetermined since the current vanishes in both configura¬ 
tions due to the triplet blockade. The parameters are the 
same as in Fig. [2l 


shown recently^ the cotunneling processes can lead to 
enhancement of TMR in this transport regime. 

Since the TMR takes well-defined values for parame¬ 
ters corresponding to plateaus in the current, it is pos¬ 
sible to find some approximate analytical formulas for 
the TMR. This can be done assuming very low tempera¬ 
tures when the Fermi functions can be replaced by step 
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functions. The formula for TMR in the transport regime 
corresponding to 5 /Ulr = and describing the TMR at 
the plateau 'URRj2<eV < ‘^Urr/2 is given by 


TMR = 


8(1+ £^) 

(7 + 9e^)l-p2’ 


( 12 ) 


while the TMR for bias voltages eV < —UrrI^. can be 
approximated by 


TMR = 


2(£^-1) 2p2 

(3e^-l)l-p2- 


(13) 


For 6/Ulr = 2, the first coefficient, 8(1 + sa)/{7 + 9^^), 
is very close to unity, while the second coefficient, 2{sa — 
1)/(3^a — 1)5 is close to zero. These two distinct values 
can be clearly seen in Fig. [6l On the other hand, for 
^ = 0 and for negative bias voltage, —UlrI2 — Ts/V^ < 
eV < —ULR/2-\-rs/V^j the TMR has a plateau of width 
V^F S' and is given by 

TMR=--^, (14) 

SI — 


while for eV < —UlrI2 — rs/\/2, the TMR reads 

1 27p 

For positive bias voltage and for (5 = 0, the TMR exhibits 
a plateau for UlrI2 — Fs/v^ < eV < UlrI2 of width 
Fs /\/2, at which it is given by 

= (16) 

Note that the TMR is always positive and smaller than 
2p^/(1 — see Fig. [6l 



eV/U,, 


FIG. 7. (Color online) The logarithm of the Andreev current 
in the parallel configuration for (5 = 0 as a function of the bias 
voltage for different Coulomb correlations F, as indicated. 
The other parameters are the same as in Fig. [2l 



eV/U,, 


FIG. 8. (Color online) The tunnel magnetoresistance as a 
function of the bias voltage calculated for (5 = 0 and for dif¬ 
ferent values of Coulomb correlations 17, as indicated. The 
other parameters are the same as in Fig. [2l 


C. The influence of intradot correlations 

Results presented in previous sections were obtained 
in the limit of infinite Coulomb correlations in the dots, 
so that the Andreev current was mediated only by CAR 
processes. We now relax this condition and allow for 
finite intradot Coulomb correlations, and study their in¬ 
fluence on the Andreev current and the TMR, focusing 
on the triplet blockade regime. Finite Coulomb correla¬ 
tions allow for nonzero current due to DAR processes, 
which can lead to a nonzero leakage current in the triplet 
blockade. Thus, once the current is finite, one can ana¬ 
lyze the behavior of the TMR, which is now well defined 
in the whole range of considered bias voltage. 

Before proceeding with the discussion of the TMR, in 
Fig. Owe first study the bias dependence of the Andreev 
current for zero detuning ^ = 0. To elucidate the role 
of finite intradot correlations the current is plotted in 
the logarithmic scale. This figure clearly shows how fi¬ 
nite Coulomb correlations affect the current in the triplet 


blockade regime. Since the dependence of the current is 
qualitatively similar in both magnetic configurations, we 
consider only the case of parallel alignment. One can see 
that in the (unphysical) case of infinite correlations, the 
current for eV > UlrI2 is suppressed in an exponen¬ 
tial way. Is oc exp(—eR/T). However, even relatively 
large values of U lead to finite current in the triplet 
blockade. For semiconductor double quantum dots the 
interdot correlations are typically an order of magni¬ 
tude smaller than the intradot correlations ^22 Although 
for recently-implemented Cooper pair splitters based on 
nanowire DQDs the capacitive coupling between the dots 
is even smallerit can still play a role. As can be seen 
in Fig. [71 the influence of direct Andreev reflection on the 
triplet blockade is clearly nontrivial, for experimentally 
relevant parameters it leads to relatively large leakage 
current in the triplet blockade, see e.g. the curves for 
U/Ulr < 10. 

The dependence of the TMR on the bias voltage is pre¬ 
sented in Fig. [8] for 5 = 0 and for different values of in- 
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tradot Coulomb correlations. First of all, one can see that 
the behavior of TMR for eV < UlrI‘^ only very weakly 
depends on U. Since the values of the TMR in this trans¬ 
port regime were discussed in the previous section, let us 
only focus on the range of bias voltages, eV > UlrI‘2, 
where the DQD is occupied by the triplet state. The 
TMR becomes then greatly enhanced and reaches values 
exceeding 2p^/(1 — p^). Such large values of TMR indi¬ 
cate that for finite intradot correlations not only the rate 
of DAR processes becomes considarable, but also that of 
CAR processes increases. This can be simply understood 
by realizing that with lowering U the occupation of the 
triplet state decreases at the cost of other states of the 
DQD, so that finite current due to both types of An¬ 
dreev reflection processes can flow. Since DAR processes 
are not sensitive to a change of magnetic configuration of 
the device (the two electrons tunnel always to the same, 
either left or right, lead), the TMR provides an indirect 
information about CAR processes in the system, the rate 
of which is clearly dependent on magnetic configuration. 

IV. SPLITTING ANDREEV BOUND STATES 

We now study how the transport properties of the 
DQD Cooper pair splitters change when one allows for 
finite splitting of Andreev bound states. Such splitting 
can be induced in various ways, e.g. by detuning the 
DQD levels, allowing for hopping between the dots24 or 
applying finite magnetic held. We again focus on the 
same parameter space as in previous sections, i.e. on 
the transport regime where mainly CAR processes are 
present. To gain a deeper insight of how the ABS be¬ 
come affected, let us consider only 9 states of the DQD 
(the limit of U oo). The analytical formula for An¬ 
dreev bound states’ energies can be then written as 

+ I + (A£)2 + 5^, (17) 

where = Sr — £r denotes the detuning of DQD lev¬ 
els, 7 ,^ = ±, and is given by Eq. (pTj) . Note 

that the splitting of ABS results only from the split¬ 
ting of single-electron states, while the states |+) and 
|—) are not affected. Moreover, one can also notice that 
finite level detuning can have a similar effect as fi¬ 
nite hopping between the dots when A^ = 2t. In the 
presence of either level detuning or hopping each of the 
four ABS states splits into two, which results in eight 
Andreev bound states. When external magnetic held is 
additionally present, these states split again and there 
are sixteen ABS states. One can thus expect that the 
consequences of ABS splitting will reveal as nontrivial 
features in transport characteristics. In the following we 
study the spin-resolved transport properties for finite de¬ 
tuning, hopping and. Anally, for finite magnetic field. 

Although we focus on transport regime where doubly 
occupied states play negligible role, in calculations we 
assume large but finite intradot correlations to be able 



FIG. 9. (Color online) The bias voltage and detuning depen¬ 
dence of (a) absolute value of the Andreev current and (b) 
the respective differential conductance in the parallel mag¬ 
netic configuration as well as (c) the TMR calculated for 
As/Ulr = 0.4. The parameters are the same as in Fig. [2] 
with UjllLR = 10. 


to determine the detuning and bias voltage dependence 
of the TMR in the considered parameter space. We thus 
assume U/Urr = 10, if not stated otherwise. 
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FIG. 10. (Color online) The bias voltage and detuning depen¬ 
dence of (a) absolute value of the Andreev current and (b) 
the respective differential conductance in the parallel mag¬ 
netic configuration as well as (c) the TMR calculated in the 
presence of finite hopping between the dots tjULR — 0.2. The 
other parameters are the same as in Fig. with UjllLR = 10. 


absent if the splitting of ABS is caused only by detuning 
of the DQD levels. The Andreev current, related differ¬ 
ential conductance in the parallel configuration and the 
TMR are shown in Fig. [9] in the case of finite and in 
Fig. [To] in the case of t = /S.e/2. The fact that t = /\e/2 
guarantees that the ABS excitation patterns occur for 
similar bias voltages and detunings in both cases. 

First of all, one can see that the number of peaks in 
differential conductance has doubled due to the two-fold 
splitting of ABS states. Because of that, the current as a 
function of the bias voltage for given detuning 5 exhibits 
more Coulomb steps compared to the case in the absence 
of ABS splitting. The region of the triplet blockade can 
be clearly visible for both finite and t, see Figs. |9]^a) 
and fTOlf a). The main difference with the case shown in 
Fig. [3] (absence of splitting) is the shift of the triplet 
line in the (eV, (5)-plane by a factor of the induced ABS 
splitting towards larger bias voltages. Moreover, the neg¬ 
ative differential conductance for 6/Ulr ^ 1 associated 
with spin accumulation in the doublet states is now also 
split, and one finds additional regions of negative G;^, 
see Figs. |9Kb) and fTOlfc). Note that for ^ > 1 , in the 
case of finite AIUlrS^ there are four regions of negative 
differential conductance, while for finite hopping t, there 
are only three. Similar asymmetry can be observed for 
negative detuning 6/Ulr ^ —1, where for A^ ^ 0 , one 
finds two regions of current suppression, which occur for 
eV < 0, while for t 7 ^ 0 there is only one, cf. Figs. [9] 
and uni However, the effect of negative differential con¬ 
ductance, which is associated with spin accumulation in 
doublet states, is not that spectacular as in the case of the 
triplet blockade, where the current suppression is much 
more pronounced (note the nonlinear color scale used in 
Figs. 13 and HOj). 

Although there are small differences in the behavior of 
the current and differential conductance in the case of 
nonzero A 5 and t, the behavior of the TMR is essential 
the same in both cases, see Figs. IlKc) andfTOlfc). One can 
observe a large TMR in the triplet blockade, TMR > 1, 
while in the case when doublet states are relevant for 
transport the TMR is close to TMR = 2p^/(l — p^), 
cf. Fig. [6] and Eq. ([12]). Now, however, this region is 
wider when changing the bias voltage and the splitting 
of ABS reveals as stripes in TMR as a function of the bias 
voltage and detuning, which are most visible for (5 > 0, 
see Figs.lHKc) andfTOlfc). 


B. Finite magnetic field 


A. Finite level detuning or hopping 

From Eq. (Ha, one could simply expect finite t and A 5 
to have the same effect on transport properties. This is 
however not entirely true, as we show in the following. 
The main difference results form the bonding and anti¬ 
bonding states that form in the case of finite t and are 


Here, we release the assumption about the smallness 
of external magnetic field Bg needed to switch magnetic 
configuration of the ferromagnetic leads. Now, we as¬ 
sume that this field is strong enough to induce splitting 
of the ABS states, which can be achieved using ferromag- 
nets with sufficiently large coercive field. The splitting 
caused by external magnetic field has a larger influence 
on spin-resolved transport compared to the splitting due 














11 



- 2-1012 


FIG. 11. (Color online) The bias voltage and detuning depen¬ 
dence of (a) absolute value of the Andreev current and (b) 
the respective differential conductance in the parallel mag¬ 
netic configuration as well as (c) the TMR calculated in the 
presence of external magnetic field BzIUlr = 0.2. The other 
parameters are the same as in Fig. [2] with UjULR = 10. 


to either finite or t. The Andreev current, differ¬ 
ential conductance in the parallel configuration and the 
TMR for finite are shown in Fig. [TTJ One can see 
that the triplet blockade region is rather not affected by 
magnetic field, see Fig. [TTl a). This is due to the fact 
that while finite Bz splits the components of the triplet, 
the total occupancy of all the triplet components does 


not change. However, the finite magnetic field changes 
the range of bias voltage for which the triplet blockade 
occurs. The negative differential conductance associated 
with the triplet blockade is thus also clearly visible, see 
Fig- fliTb). Interestingly, the negative differential con¬ 
ductance due to the spin accumulation in the doublet 
states in now present only for negative values of detun¬ 
ing 6/Ulr ^ —1, cf. Figs. [3] and [TTJ This can be un¬ 
derstood by realizing that while for 6/Ulr ^ —1 and 
eV < 0, splitting caused by finite magnetic field addi¬ 
tionally enhances the spin accumulation, in the case of 
^/Ulr ^ 1 and eR > 0, magnetic field diminishes the 
spin accumulation. Consequently, the current suppres¬ 
sion is reduced in the latter case, while in the former one 
it is enhanced. 

The bias voltage and detuning dependence of the TMR 
is shown in Fig. fTH c). First of all, one can notice 
a strong asymmetry with respect to the bias reversal, 
which is most visible for small eV and SIUrr ^ 1- For 
5/Ulr ^ 1 and for low bias, the DQD is empty and 
the TMR is positive in this transport regime, and rather 
symmetric around eV = 0. On the other hand, when 
\^/Ulr\ ^ 1, the DQD is singly occupied and the TMR 
reveals then a strong asymmetry with respect to the bias 
reversal. For positive bias voltage, there is a large pos¬ 
itive TMR, while for negative bias, the TMR becomes 
negative. Such asymmetry is associated with the split¬ 
ting of the doublet ground state of the DQD. It strongly 
affects CAR processes in the antiparallel configuration, 
since the Cooper pair electrons tunnel then to either ma¬ 
jority or minority spin bands of the ferromagnets. Be¬ 
cause for one bias polarization the electron occupying 
the DQD is the majority-band electron, while for the op¬ 
posite bias polarization, this electron belongs to the spin 
minority channel, it effectively leads to large asymme¬ 
try of the flowing current with respect to the bias re¬ 
versal, which is most visible in the TMR. The effect of 
sign change of the TMR can be even more pronounced 
in the case of 6/Ulr ^ —1, where the DQD is occupied 
by two electrons. For eV > 0, one then finds TMR ^ 2, 
while for eV < 0, TMR ^ —1/2. The mechanism lead¬ 
ing to this asymmetry is similar to that described above. 
Note, however, that in the blockade regions the current, 
and thus the TMR, can be still modified by cotunneling 
processes 

Out of the Coulomb blockade region, the changes in 
TMR are not that spectacular, however, there are still 
considerable differences compared to the case of = 0. 
First of all, negative values of the TMR can be observed 
for eV < —Ulr and 6/Ulr ^ —1- Moreover, although 
in the triplet blockade regime, one observes a large pos¬ 
itive TMR, similarly as in the absence of magnetic field, 
cf. Figs. [6] and Fig. [TTj^c), for finite Bz there is an addi¬ 
tional region of an enhanced TMR. It occurs for voltages 
around eV ^ {6 ^ Ulr)/‘^ and is associated with the 
splitting of the triplet components. Such splitting in¬ 
creases the difference in the currents in both magnetic 
configurations, yielding TMR > 1. 
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5/Ulr 

FIG. 12. (Color online) The absolute value of the Andreev 
current (a) calculated for the parallel (Ig) magnetic config¬ 
uration, the corresponding differential conductance (b) and 
the tunnel magnetoresistance (c) as a function of detuning 
S = 2s Ulr and the applied bias eV. The parameters are 
the same as in Fig. [2] with UjULR = 2. 


V. RESULTS IN THE FULL PARAMETER 
SPACE 

To complete the analysis of Andreev transport through 
DQD Cooper pair splitters, in this section we extend the 
discussion to the full parameter space. Figure [12] presents 
the bias voltage and detuning dependence of the abso¬ 


lute value of the Andreev current and the corresponding 
differential conductance in the parallel magnetic configu¬ 
ration as well as the TMR. Transport characteristics are 
now symmetric with respect to the particle-hole symme¬ 
try point of the DQD Hamiltonian, Sph = —Ulr — U /2 
{^ph = —Ulr — U)^ with immediate sign change of the 
bias voltage, Is{eV,S > Sph) = -Is{-eV,S < Sph), 
Gs{eV,6 > 6ph) = Gs{-eV,6 < 5ph), and TMR(eU,^ > 
5ph) = TMR(—eU, 5 < 5ph), see Fig. [121 Note that to en¬ 
able direct comparison with previous results we still plot 
transport characteristics as a function of S = 2s -\-Urr- 
Moreover, we assumed relatively large capacitive cou¬ 
pling between the two dots, U = 2UlrU^^^ able to show 
transport properties in the whole range of detuning S in 
a single panel and not to obscure the features discussed 
previously, which occur around (5^0. However, results 
are qualitatively the same for larger intradot correlations, 
as checked numerically (not shown), the main difference 
is in the distance between the resonances occurring now 
for 6 /Ulr ~ —1 and S/Ulr ~ —5, see Fig. [121 which 
increases with increasing U. 

Due to finite Coulomb correlations in the dots, there 
are more Andreev states available for transport, which 
generally reveals as steps in the bias dependence of the 
current and corresponding peaks in the differential con¬ 
ductance. One can clearly see the regime of the triplet 
blockade, which occurs for {S + Ulr)/‘^ ^ eV ^ {S -\- 
Ulr + ^^U) and ^ 0 or (5 < —2U — 2Ulr- Note that there 
is a relatively large leakage current in the triplet block¬ 
ade due to finite intradot correlations, which allow DAR 
processes to participate in transport. Interestingly, there 
are also another regions where the current becomes sup¬ 
pressed and negative differential conductance occurs, see 
Figs. m a) and (b). To present and discuss these effects 
in more detail, let us show the relevant cross-sections of 
Fig. dll Since all transport features display an appropri¬ 
ate symmetry with respect to ^ = 6ph, in the following 
we will only analyze the results for 6 > 5ph. 

The current and differential conductance in both mag¬ 
netic configurations as well as the resulting TMR in the 
case of 6/Ulr = —3 are shown in Fig. [TS] First of all, 
we note that transport characteristics are now symmetric 
with respect to the bias reversal. This is due to the fact 
that the presented data were obtained for the particle- 
hole symmetry point of the model S = SpR. Moreover, 
one can see that the current does not increase in a mono¬ 
tonic way, see Fig.[T3]^a). For \eV\/ULR ^ 1, the Andreev 
current becomes suddenly suppressed and the system ex¬ 
hibits a pronounced negative differential conductance, 
see Fig. mb), which is present in both magnetic con¬ 
figurations. The decrease of the current is related with 
an enhanced occupation of doublet states. More pre¬ 
cisely, for positive bias voltage cV/Ulr ^ 1 and in the 
antiparallel configuration, the occupation of the states 
(It, (i) —|(i, t))/v^ and (| t, d) — \d^ \))/^/2 is close to unity, 
which decreases the rate for both DAR and CAR Andreev 
processes. Similar situation also holds for negative bias 
voltage, but now the one-electron doublet states become 
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FIG. 13. (Color online) The Andreev current (a) and the 
corresponding differential conductance (b) in both magnetic 
configurations together with the tunnel magnetoresistance (c) 
as a function of applied bias eV for SIUlr = —3 ((5 = Sph)- 
The other parameters are the same as in Fig. 1121 


occupied. Consequently, the Andreev current becomes 
suppressed and the system exhibit negative differential 
conductance. However, with further increase of the bias 
voltage, \eV\/ULR the occupation of the above dou¬ 
blet states decreases and the current raises again, chang¬ 
ing then monotonically with the bias voltage. 

The above-described behavior is present in both mag¬ 
netic configurations, however, in the parallel configura¬ 
tion there is a strong spin accumulation in the doublet 
states and the current is more suppressed compared to 
the antiparallel configuration. This is reflected in the be¬ 
havior of the TMR on the bias voltage, which is shown in 
Fig. me) . At low voltage the TMR is negligible, which 
is related to the fact that the system is occupied by the 
singlet state (a(|0, 0) — jd, d)) -h/3(| t, i) — | t))? ther¬ 

mally activated transport through this state is insensitive 


FIG. 14. (Golor online) The same as in Fig. fT3l calculated for 
S/Ulr = - 2 . 


to the change in magnetic configuration. With increasing 
eR, the TMR increases at the first step in the current. In¬ 
terestingly, for voltages where the current is suppressed, 
the TMR takes large values due to spin accumulation in 
the doublet states, but then drops again with next step 
in the current when the voltage is increased further on. 

When moving away from the particle-hole symmetry 
point, 6 = 5phi there is a large change in the transport 
characteristics, see Fig. [121 First of all, a pronounced 
asymmetry with respect to the sign change of the bias 
voltage occurs. Moreover, transport characteristics be¬ 
come more complex, since the number of negative differ¬ 
ential conductance regions increases and one can also find 
transport regimes where negative TMR occurs. Here, 
let us discuss in somewhat greater detail the case of 
^/Ulr = —2, which is presented in Fig. [HI 

The Andreev current as a function of bias voltage is 
shown in Fig. [T^ a). Its bias dependence is not mono¬ 
tonic irrespective of bias polarization and magnetic con- 
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figuration of the device. For both positive and nega¬ 
tive bias voltage, there are two regions of current sup¬ 
pression accompanied with respective negative differen¬ 
tial conductance, see Fig. HU^b). For positive bias, the 
current first decreases once eVjULR ~ 3/2, which is as¬ 
sociated with enhanced occupation of the doublet state 
[| t, d) — \d^"\')]/^/2. In this transport regime DAR pro¬ 
cesses are suppressed and transport is mainly due to CAR 
processes. With increasing the bias voltage. Is increases 
for eV/U lr ~ 5/2 to drop again once eV/U lr ~ 4, where 
the DQD becomes mainly occupied by the state |d, d). 
Then, the rate of both Andreev reflection processes be¬ 
comes decreased. Similar features can be observed for 
negative bias voltage and the mechanism leading to cur¬ 
rent suppression and negative differential conductance is 
basically the same. The first negative differential conduc¬ 
tance occurs due to enhanced occupation of the doublet 
state [|4,,0) - | 04 )]/V 2 (for -1/2 > cV/Ulr > -3/2), 
while in the second suppression region the DQD is in the 
state |0,0) (for eVjULR ^ —2). 

Interestingly, in regions of current suppression due to 
enhanced doublet occupation, the TMR exhibits rather 
large values, see Fig. He). This indirectly confirms that 
the main role is played by CAR processes, which greatly 
depend on magnetic configuration of the system, contrary 
to DAR processes. Consequently, one observes a large 
positive TMR effect. On the other hand, in blockade re¬ 
gions due to enhanced occupation of either empty or fully 
occupied DQD, the Andreev current depends very weakly 
on the magnetic configuration of the system, which im¬ 
plies that CAR processes play a minor role in transport. 

One can also note that CAR processes become rele¬ 
vant not only in blockade regions, but also for 1/2 < 
eVjULR < 3/2, that is at the first plateau for positive 
bias voltage, see Fig. He). This is related with nonequi¬ 
librium spin accumulation in the parallel configuration, 
due to which the occupation of the triplet component 
|t,t) becomes enhanced and the Andreev current drops. 
This triplet blockade is absent in the antiparallel config¬ 
uration, which leads to large difference in the currents in 
both configurations and thus to large TMR effect. 

Another interesting feature is the negative TMR, 
which occurs at very low positive bias voltage, see 
Fig. [Mf cb In this transport regime the Andreev pro¬ 
cesses are suppressed due to the fact that the DQD is 
occupied by the doublet state [| 0) — |0 ,|)]/a/ 2 with 

almost unit probability. This occupation probability is 
slightly lower in the parallel configuration for positive 
bias voltage at the cost of small but finite occupation 
of the state a[\ t, 0) + |0, t)] + P[\d) + jd ,']')]. It is this 
state that allows for finite Andreev current in the parallel 
configuration, giving rise to negative TMR effect. 
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FIG. 15. (Color online) Fidelity F in the parallel magnetic 
configuration as a function of detuning 6 and the bias volt¬ 
age eV calculated for (a) As = 0 and (b) AsIUlr = 0.4. 
The dark (green) region corresponds to F = I. The other 
parameters are the same as in Fig. [2l 

Cooper pair is split and each electron tunnels to differ¬ 
ent arm of the device in a CAR process, we again focus on 
the transport regime where DAR processes are excluded 
by assuming infinite intradot Coulomb correlations. 

To analyze the fidelity let us consider the Werner 
stated which has the following form, 

Tr(F)=F|5)(5| + (l-F)Fd|^, (18) 

where F denotes Werner fidelity. For F < 1/2, the 
Werner state is unentangled, whereas for 1/2 < F < 
I, there exist purification protocols, which can extract 
states with arbitrary large entanglement. Werner fidelity 
for the considered system is given by the formula^ 


VI. ENTANGLEMENT FIDELITY F = — -—, (19) 

Fs + Ft 

In this section we study the entanglement fidelity be- where Pt = and P^ are occupation proba- 

tween split electrons forming a Cooper pair. Since the bilities for the triplet and singlet states, respectively. The 
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latter can be expressed by P± = {±\p\±) as 


Ps 




P_ 




( 20 ) 


with a = ( 1 / a / 2)^1 — S/{2sa) and /3 = 

(l/\/2)i/l + 5/{2eA)‘ Roughly speaking, when 

Ps ^ Pti fidelity reaches its maximal value F ^ 1, 
while in the opposite situation it is suppressed. 

Fidelity in the parallel magnetic configuration as a 
function of S and eV is shown in Fig. [15] (a) in the absence 
of detuning, = 0, and (b) for Ae/Ulr = 0.4. One can 
clearly observe transport regimes where F is either equal 
to one or zero. More specifically, for eV < \ {S PUlr)/‘2\^ 
one has Ps ^ Pti and fidelity reaches its maximal value 
with F ^ 1. Thus, the transmitted pairs of electrons can 
be considered as entangled. On the other hand, for bias 
voltages, eV > |((5 + FLi?)/2|, the situation is just oppo¬ 
site and one obtains F ^ 0. Consequently, the considered 
DQD setup guarantees that, by properly tuning the de¬ 
vice parameters, fully entangled pairs of electrons can be 
extracted from the superconductor and transmitted into 
normal leads. This effect is insensitive to the value of the 
leads’ spin polarization p and the magnetic configuration 
of the system (results not shown). It is also interesting to 
notice that fidelity provides information about the flow¬ 
ing current. The current does not flow due to the triplet 
blockade {Ps = 0 and Ft = 1), i.e. when F = 0, cf. 
Figs. [2] and [151 

The situation becomes more complex when finite de¬ 
tuning of DQD levels is present. In this case the map of fi¬ 
delity possesses richer structure, see Fig.[T5lb). The main 
difference is in the splitting of the line, eV ^ {SpUlr)/‘2, 
along which the fidelity varies between zero and one. 
Now, one obtains F = 1 for smaller bias voltages for 
given detuning S > —UlrI‘2^ compared to the case of 
Ae = 0. In fact, the voltage is smaller by a factor of 
level splitting A^. Moreover, with increasing the bias 
voltage and for S > —PlrI‘^^ F does not drop to zero 
immediately, but becomes suppressed in a nonmonotonic 
way in the transition region, eV ^ {6 Urr)/2 ± A^, of 
width 2 A 5 , see Fig. [T5fbb Finally, we note that similar 
splitting of the transition line separating the regions with 
F = 0 and F = 1 also occurs in the case of finite hopping 
between the dots or finite magnetic field. 


coupled to a common 5-wave superconductor and each 
dot coupled to its own ferromagnetic lead. The calcula¬ 
tions were performed with the aid of the real-time dia¬ 
grammatic technique, assuming weak coupling between 
DQD and ferromagnets and taking into account the se¬ 
quential tunneling processes. We have analyzed the bias 
voltage and DQD level dependence of the Andreev cur¬ 
rent and the differential conductance in the parallel and 
antiparallel configurations, as well as the resulting tunnel 
magnetoresistance. 

In the case of infinite correlations in the dots, we have 
discussed the behavior of spin-dependent characteristics 
in the transport regime where only crossed Andreev re¬ 
flection processes are possible. For certain DQD levels’ 
configuration and applied bias voltage, the current is then 
suppressed due to the triplet blockade. We showed that 
even relatively large intradot correlations can lead to fi¬ 
nite leakage current in the triplet blockade. We found an 
enhanced TMR in the triplet blockade, which indicates 
the role of CAR processes in transport. We have also 
analyzed the effect of splitting the Andreev bound states 
by either finite DQD level detuning, finite hopping be¬ 
tween the dots or finite magnetic field. While in the first 
two cases each Andreev bound state becomes split into 
two, finite magnetic field further splits the ABS, resulting 
in more complex transport characteristics, with negative 
differential conductance and negative TMR occurring in 
certain transport regimes. 

Moreover, assuming finite correlations in the dots, we 
have studied transport properties in the full parameter 
space, where both DAR and CAR processes are relevant. 
We found transport regimes where additional current 
suppression accompanied with negative differential con¬ 
ductance occurs. These suppression regimes are due to 
enhanced occupation of certain many-body DQD states, 
which diminishes the rate of either CAR or DAR pro¬ 
cesses, depending on transport region. 

Finally, in the CAR transport regime we have also ana¬ 
lyzed the entanglement fidelity between electrons forming 
Cooper pairs. We showed that the fidelity of split Cooper 
pair electrons can be tuned by bias and gate voltages and 
for certain parameters F can reach unity. Consequently, 
DQD-based Cooper pair splitters, by properly tuning the 
device parameters, can be sources of fully entangled pairs 
of electrons that are extracted from superconductor and 
transmitted to normal leads. 


VII. CONCLUSIONS 

In this paper we have studied the spin-resolved An¬ 
dreev transport through double quantum dot-based 
Cooper pair splitters with ferromagnetic leads. The con¬ 
sidered device consisted of two single-level quantum dots 
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